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Abstract: This study aimed to identify the assemblages (or subassemblages) of Giardia duodenalis by using normal or 
nested PGR based on 4 genetic loci: glutamate dehydrogenase (gdh), triose phosphate isomerase (tpi), p-giardin (bg), 
and small subunit ribosomal DNA (18S rRNA) genes. For this work, a total of 216 dogs' fecal samples were collected in 
Guangdong, Ghina. The phylogenetic trees were constructed with MEGA5.2 by using the neighbor-joining method. Re- 
sults showed that 9.7% (21/216) samples were found to be positive; moreover, 10 samples were single infection (7 iso- 
lates assemblage A, 2 isolates assemblage G, and 1 isolate assemblage D) and 1 1 samples were mixed infections where 
assemblage A was predominant, which was potentially zoonotic. These findings showed that most of the dogs in Guang- 
dong were infected or mixed -infected with assemblage A, and multi-locus sequence typing could be the best selection 
for the genotype analysis of dog-derived Giardia isolates. 
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Giardia duodenalis is one of the most common intestinal 
protozoan parasite, which is ubiquitous in mammals [1]. In 
Asia, Africa, and Latin America, approximately 200 million 
people present symptomatic giardiasis and 500,000 new cases 
are diagnosed per year [2]. 

G. duodenalis is considered as a species complex, comprising 
at least 8 distinct genetic groups, referred to as assemblages A 
to H [3]. To date, only assemblages A (subtypes I and II) and B 
(subtypes I and IV) have been associated with human infec- 
tions, but are also found in a number of other mammalian 
hosts [4]. Genetic assemblages C, D, E, F, G, and H appear to 
be host-restricted to domestic animals and wild animals [5]; 
furthermore, assemblages C and D infect dogs and wild cami- 
vores, assemblage E infects cattle and other hoofed livestock, 
assemblage F infects cats, assemblage G infects rats [6], and as- 
semblage H infects marine mammals [7]. 

In China, there are approximately 28.5 million cases per 
year The majority of giardiasis cases remain unreported [8]. 
However, Yong et al. [9] have obtained 8 G. duodenalis isolates 
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from Anhui Province, China (all from purified cysts); 4 iso- 
lates were grouped into assemblage A and 4 into assemblage B. 
In Henan Province, China, 18 Giardia-positive specimens were 
characterized by sequence analysis of the triose phosphate 
isomerase (tpi) gene into assemblages A (8 belonging to AI 
and 4 All) and B (belonging to 6 new subtypes) found in 12 
and 6 specimens, respectively [10]. What's more, 5 of 23 iso- 
lates from dogs were characterized as assemblage AI and 18 as 
assemblage D [11]. Dogs have been considered as the main 
potential reservoirs because of their close contacts with hu- 
mans, but direct evidence for animal-to-human transmission 
is still scarce [1,12]. So far, most studies have been mainly 
based on analysis of the small subunit ribosomal locus (18S 
rRNA), |3-giardin (bg), glutamate dehydrogenase (gdh), elon- 
gation factor 1-alpha (ef-1), and triose phosphate isomerase 
(tpi). Only bg, gdh, and tpi genes were sufficiendy discrimina- 
tory to type at the level of subassemblage, 18S rRNA and ef-1 
permit identification only at the level of the assemblage, which 
are much more conserved than other three loci [6,13]. 

To date, due to the lack of high-resolution genotype identifi- 
cation based on 1 gene, a multi-locus sequence-typing scheme 
is needed on the molecular identification of zoonotic assem- 
blages in dogs. Therefore, the aim of the present study was to 
genetically characterize isolates of G. duodenalis from dogs in 
Guangdong, China using analysis of tpi, gdh, bg, and 18S 
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rRNA loci combined. 

A total of 216 dog fecal samples were obtained from 4 dif- 
ferent shelters in Guangdong, China during June 2012 to Oc- 
tober 2013. Two shelters were located in Shunde (urban) and 
Sanshui (suburban) districts of Foshan City, Guangdong and 
2 in Baiyun (urban) and Conghua (suburban) districts of 
Guangzhou City, Guangdong. Thirty fecal samples were col- 
lected in Shunde distria, 55 in Sanshui district, 80 in Conghua 
district, and 51 in Baiyun district. All the dogs were abandoned 
or stray dogs and had stayed in the shelters for at least 1 
month. Fresh fecal samples from each dog were transported 
back to the Parasitology and Parasitic Diseases Laboratory 
College of Veterinary Medicine, South China Agricultural Uni- 
versity (SCAU) on the same day of collection, preserved in 
2.5% potassium dichromate and kept at 4°C. Direct micro- 
scopic examination was done after centriftigation of 1.5 g feces 
by the flotation technique with saturated glucose 

Positive samples examined by microscopy were pretreated 
with 5 cycles of hearing at 100°C for 5 min and immediate 
freezing at -80°C for 5 min. DNA was extracted directly from 
fecal samples using a commercial DNA extraction kit (QIAamp 
DNA Stool Mini Kit, QIAgen, Hilden, Germany) according to 
the manufacturer's instmcdons. DNA samples were stored at 
-20°C until use. 

The IBS rRNA, bg, and gdh genes were amplified individu- 
ally according to the protocols described previously [6,14,15]. 
The tpi primers, tpil (5' AACGCAATCACrGTATCT 3') and tpi2 
(5' CAATGACAACCrCCTTCC 3') were used for the primary 
amplification, and tpi3 (5' CTTCATCGGCGGTAACn 3') and 
tpi4 (5' GGCACGCnAGCCrrcrr 3') were used for the nested 
amplification. Identical conditions for tiie primary and sec- 
ondary amplification of the tpi gene fragment were as follow- 
ing: 35 cycles (94°C for 30 sec, 45°C for 30 sec, and 72°C for 45 
sec) in a T-personal thermocyder, witii an initial preheating at 
94°C for 5 min and a final extension at 72°C for 7 min. PGR 
products were analyzed after electtophoresis in 1.5% agarose 
gels and stained with 0.2 g/ml of ethidium bromide and visu- 
alized on a UV transilluminator Positive amplicons were puri- 
fied and sequenced at an ABI 377 automated DNA sequencer 
(BigDye Terminator Chemistry). Obtained sequences were 
aligned by Clustal X (2.21) using default parameters, and ana- 
lyzed by MEGA5.2 software. Obtained nucleotide sequences 
have been submitted to the GenBank™, under accession num- 
bers (KJ027399 to KJ027403 and KJ027407) for 18S rRNA; 
(KJ027417 to KJ027423, KJ027411, KJ027413, KJ0274134, 



KJ027415, and KJ027425) for bg, (KJ027426 to KJ027433, 
and KJ027438) for gdh, and (KJ027448 to KJ027463) for tpi. 

A 334 bp fragment of the tpi gene was obtained from 16 
dog's isolates; the sequence analysis revealed assemblage Al in 
15 isolates and assemblage D in 1 isolate. The amplification of 
a 515 bp fragment of the bg gene was obtained from 12 dog's 
isolates; sequence analysis revealed assemblage Al in 1 isolate, 
assemblage C in 3 isolates, and assemblage D in 8 isolates. 
The amplification of a 530 bp fragment of the gdh gene was 
obtained from 8 dog isolates; sequence analysis revealed as- 
semblage Al in 1 isolate, assemblage C in 3 isolates, and as- 
semblage D in 4 isolates. The amplification of a 292 or 291 bp 
fragment of the 18S rRNA gene was obtained from 6 G. duode- 
nalis isolates, the sequence analysis revealed assemblage A in 4 
isolates, assemblage C in 1 isolate, and assemblage D in 1 iso- 
late (Table 1). 

A single infection with assemblage A was found in 7 sam- 
ples, where 1 was typed by 3 loci, 3 were typed by 2 loci, and 3 
was typed by 1 locus. A single infection with assemblage C was 
found in 2 samples, where 1 was typed by 2 loci and 1 was 
typed by 1 locus. A single infection with assemblage D was 
found in 1 sample that was typed at 3 loci. Among mixed in- 
fections, 2 assemblages were detected in 11 samples (A and D 



Table 1 . Genotyping data of Giardia intestinalis isolates from dogs 
in Guangdong, China at tlie 4 loci tpi, gdli, bg, and 18S rRNA 
Neg, negative 



Isolates 


tpi 


gdh 


bg 


18S rRNA 


DBY14 


Al 


Neg 


Neg 


A 


DBY23 


Al 


D 


D 


Neg 


DBY33 


Al 


Neg 


D 


Neg 


DBY37 


Al 


C 


Neg 


Neg 


DBY42 


Al 


Neg 


D 


Neg 


DBY50 


Al 


Neg 


Neg 


A 


DCH08 


Neg 


C 


C 


Neg 


DCH25 


Neg 


Neg 


C 


Neg 


DCH62 


Al 


D 


D 


Neg 


DCH65 


Al 


Neg 


D 


Neg 


DCH67 


Al 


C 


Neg 


C 


DCH71 


Al 


D 


Neg 


Neg 


DCH78 


Al 


Neg 


Neg 


Neg 


DSS10 


Al 


Neg 


Neg 


Neg 


DSS17 


Al 


Neg 


D 


Neg 


DSS18 


Neg 


Neg 


Al 


A 


DSS30 


Al 


Neg 


Neg 


Neg 


DSS48 


Al 


Al 


Neg 


A 


DSD05 


D 


Neg 


C 


Neg 


DSD12 


Neg 


D 


D 


D 


DSD28 


Neg 


C 


D 


Neg 
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in 7, D and C in 2, A and C in 2), more than 2 assemblages in 
the same sample were not found (Table 2). 

The phylogenetic analysis was performed by using the 
neighbor-joining method based on 4 gene loci. Result showed 
that assemblages A and D were significantly separated at the 
tpi gene. Assemblage A subgroups remained separated in the 
tpi gene. All sequences of assemblage A were placed in 1 clus- 
ter with AI sequence reference (GU5 64275) separating from 
All and AIII in the neighbor joining tree (Fig. lA). Assemblag- 
es A, C, and D were significantly separated at the gdh gene, 
and subgrouping was observed among assemblages A in the 



gdh analysis. The neighbor-joining tree placed the representa- 
tive sequence (DSS48) in 1 cluster with AI sequence reference 
(EF685696) with high bootstrap support (Fig.lB). Assem- 
blages A, C, and D were significantly separated at the bg (Fig. 
IC) and 18S rRNA genes (Fig. ID). At bg gene, the sequence 
(DSS18) were placed in 1 cluster with AI sequence reference 
(AB469365) separating from All and AIII. Phylogenetic Uee of 
the tpi, gdh, bg, and 18S rRNA genes by the maximum-likeli- 
hood method were similar to that by the neighbor-joining 
method. 

The homologous analysis was performed by using nudeo- 



Table 2. Assemblages and subassemblages of 21 positive samples from dogs at 1 , 2, or 3 loci 



Locus/Loci 






Assemblages and sub-assemblages 






1 locus 


tpi 
Al(3) 


gdh 


bg 
C(1) 


18S 






2 loci 


tpi, bg 


tpi, gdii 


tpi, 18S 


bg, gdh 


bg, 18S 


gdh, 18S 




AI, D (4) 


AI, C(1) 


AI, A (2) 


C,C(1) 


AI, A(1) 






D,C(1) 


AI, D(1) 




D,C(1) 






3 loci 


tpi, bg, gdii 


tpi, gdh, 18S 


bg, gdh, 18S 


tpi, bg, 18S 








AI, D, D (2) 


AI,AI, A(1) 
AI, C, C{1) 


D, D, D(1) 









Results in "1 locus" labeled are from samples which were amplified only at this locus; Results in 2 and 3 loci are for samples with assemblages and 
subassemblage information at more than 1 locus. For example, in the "2 loci" labelled tpi and bg first sample is AI,D (4), which indicates this sample 
had subassemblage AI for tpi and assemblage D for bg, and the number 4 in parentheses indicates that there were 4 samples with this result. 



Table 3. Summary of nucleotide variations of Giardia intestinalis A genotype subassemblages in dog isolates at the bg, tpi, and gdh genes 



Alignment position 



Reference sequences 


Bg 

AI-AB469365 
AII-AY072723 
AIII-AY072724 


68 
T 
C 
G 


317 
T 
C 
G 


327 
G 

T 


335 
T 

G 


473 
G 
T 
T 






















Dog 


DSS18 


G 


G 




























Reference sequences 


Gdh 


63 


118 


162 


180 


258 


312 


366 


390 


420 


426 


429 


453 


461 








AI-EF685696 


G 


T 


T 


G 


T 


C 


G 


G 


T 


T 


T 


C 


G 








AII-EU362962 


A 




G 


T 


G 


T 


T 


T 


C 


G 


G 


T 


T 






Dog 


DSS48 
































Reference sequences 


Tpi 


25 


50 


58 


73 


76 


82 


85 


94 


98 


109 


127 


139 


146 


154 


196 246 




AI-GU564275 


T 


G 


T 


G 


G 


C 


G 


T 


G 


G 


A 


A 


A 


A 


C T 




AII-GU564277 












T 


T 




T 


A 


G 


G 




G 


T 




AIII-EU781002 






G 


T 


A 






G 
















Dog 


DSS10 
DBY50 
DBY23 

DSS30/DSS17/ 
DGH78/DCH71/ 
DGH67/ DCH65/ 
DGH62/ DBY42/ 
DBY37/ DBY33/ 
DBY14/DSS48 




A 






















G 




G 



Nucleotide substitutions are numbered from the start of the gene which was amplified. Dots indicate identity to the AI reference sequence. 
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72 



95 



83 



65' 



DBY50 

♦ DBY42 

♦ DCH71 

A1-GU564275 

♦ DCH78 
-♦ DSS10 
♦dCH62 

36 ♦ DSS30 

♦ DCH67 
-♦ DBY23 
♦dSS17 

94 ♦dCH66 
DBY37 
'DBY33 
' DBY14 
•DSS48 
A2-GU564278 
62 'a2- GU564277 
A3-EU781002 
A3-DQ660648 

F-JF968100 

E-JN 162347 

D-DQ220289 
DSD05 
B-FJ560571 



100 



0.01 



44 

58 A 
100 



♦ DSD12 

♦ DSS17 

♦ DCH65 

♦ DCH62 

♦ DBY42 

♦ DBY23 

♦ DBY33 

♦ DSD28 
D-AY546647 



86 



82 



DSD05 
gg I C-AY545646 
48lr* DCH08 

40^ ♦ DCH25 

B-AY072725 



79 



■ F-AY647264 



51 



100 



35 J* DSS18 
A1-X86968 
A3-AY072724 



54 A2-AY072723 



0.01 



B 



92 
100 



79 



♦ DBY23 

♦ DSD12 

♦ DCH62 

♦ DCH71 

— D-EF507612 



100 



♦ DSD28 

a3| C-EF507621 

'♦ DBY37 
63 M ♦ DCH67 
69 ♦ DCH08 
— B-EF685682 



99 



E-AY1 78740 



54 



100 



97 



F-AB569381 

AII-EU362962 
AII-EU362957 
AI-EF685696 



99 ♦ DSS48 



I 1 

0.01 



D 



♦ DBY50 



6 
7 
17 



■ B-AB569368 



♦ DSS18 

♦ DBY14 

♦ DSS48 
A-AB469363 



- F-DQ836339 

- E-AB569370 



89 



82 



DSD12 

' D-AB569371 
DCH67 



651 C-AB569372 



I 1 

0.002 



Fig. 1 . Phylogenetic relationships of Giardia duodenalis assemblages from dogs inferred by the neighbor-joining analysis of 4 gene loci. 
(A) Triose phosphate isomerase (tpi) nucleotide sequences. (B) Glutamate dehydrogenase (gdh) nucleotide sequences. (C) P-giardin (bg) 
nucleotide sequences. (D) Small subunit ribosomal DNA (1 8S rRNA) nucleotide sequences. Bootstrap values obtained from 1 ,000 repli- 
cates are indicated on branches in percentage. 



tide sequence based on tpi, gdh, and bg gene loci. At tpi gene, 
all sequences of assemblage A showed complete homology to 
AI except for DSSIO (G-A at position 50), DBY50 (A-G at posi- 
tion 146) and DBY23 (T-C at position 246) of dog isolates. At 
gdh gene, DSS48 sequence showed complete homology to 
the AI reference sequence (EF685696). At bg gene, DSS18 se- 
quence showed complete homology to AI reference sequence 



(AB469365) (Table 3). 

The current study present the first genotyping by multi-locus 
sequence typing of G. duodenalis isolated from dogs in Guang- 
dong, China. Most of the dogs in this study were infected or 
mixed-infected with assemblage A, which is different from 
most previous studies. Recently Feng and Xiao [8] have stated 
that 1,049 of 1,563 (67%) G. duodenalis isolates from dogs 
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were either assemblage C or D in the world. This phenome- 
non may be related to the living environments of the dogs, as 
they lived in shelters, the first infection by assemblage A may 
predominate over the shelters latterly [ 1 3] . On the other hand, 
dogs from shelters which aowded breeding show higher rate 
of dog-to-dog transmission and so do the young animals than 
older ones [16]. Epe et al. [17] have reported an average rate of 
24.8% in dogs through Europe, which rises to 34.7% in shel- 
ter dogs. 

In our study, the tpi locus show the highest amplification 
rate in dogs, where most of samples were typed as assemblage 
Al except 1 (DSD05) which was typed as assemblage D. The 
sensitivity of tpi gene seemed to be high in amplifying assem- 
blage A which has zoonotic potential. Two loci gdh and bg ap- 
peared more likely to amplify assemblage C and D in dogs, 
which is in coherence with Scorza et al. [18]. Moreover, the 
gdh and bg gene also amplified 1 assemblage Al in dogs re- 
speaively, which have zoonotic potential with certainly. Our 
results showed the presence of potentially zoonotic assemblag- 
es in dogs in Guangdong, China, which is consistent with Li et 
al. [11]. Comparing the conserved marker (18S) with variant 
marker (bg, gdh, and tpi) in genotyping and the amplification 
rate, the PCR amplification at the 18S rRNA loci had the poor- 
est performance, but with longer fragment (about 292 bp) in 
this study than most of tiie previous studies. To date, almost 
half of the studies performed were based on the analysis of a 
short fi^agment (130 bp) of 18S rRNA locus [19]. Furtiiermore, 
studies on G. duodenalis have carried out for several years in 
our lab, and no assemblage B was found in dog's isolates from 
Guangdong, South China. The assemblage B isolates might 
not exist in this region. Thus, more researches are stiongly re- 
quired to clarify the presence/absence of assemblage B. 

The high frequency of multiple/mixed infections found in 
this study has been inaeasingly reported in other countries on 
multi-locus studies in both humans and mammalians includ- 
ing dog, cat, cattie, goat, sheep, pig, and wildlife specimens 
[6,8,12,13,20-22[. Eleven samples were mixed infections 
where assemblage A was predominant in die study. It may be 
atuibuted to mixed infections or allelic sequence heterozygos- 
ity (ASH). Mixed infections can happen when a host ingests 
Giardia cysts of different genetic profiles or subsequent infec- 
tion of an infected host by genetically different Giardia cysts. 
This is especially common in areas where giardiasis is endemic 
[3,6,12,20,23[, but ASH can also account for this finding and 
sexual reproduction may influence ASH levels [12]. However, 



the contribution between mixed infections and ASH in the 
complexity of genotyping is currentiy uncertain. 

In conclusion, assemblage A was most frequent among in- 
fected dogs, which is potentially zoonotic genotype. Moreover, 
dog giardiasis seems to be more likely dog-to-dog transmis- 
sion, but humans may be involved in a zoonotic cycle of 
transmission, where human-to-human transmission is much 
more important than animal-to-human transmission. As well, 
more studies on the contributing factors between mixed infec- 
tions and ASH should be seriously taken place. The choice of 
gene loci can influence the ampliflcation results, hence, detec- 
tion methods targeting loci can be extremely useful when a 
common source of contamination is certainly involved. Thus, 
multi-locus sequence typing could be considered for the geno- 
type analysis of dog-derived Giardia isolates. 
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